Kondo effect in double quantum dots with inter-dot repulsion 
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We investigate a symmetrical double quantum dot system serially attached to the leads. The emphasis is on 
the numerical analysis of finite inter-dot tunneling in the presence of inter-dot repulsive capacitive coupling. 
The results reveal the competition between extended Kondo phases and local singlet phases in spin and charge 
degrees of freedom. The corresponding phase diagram is determined quantitatively. 
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Quantum dots: 1 ' 2 - 3 - 4 provide in addition to various applica- 
tions in proposed spintronics 5 and quantum computation 6 - 7 
devices a playground for studying phenomena known in bulk 
condensed-matter systems. Specifically, the Kondo effect was 
found to play an important role in single 8 and double quan- 
tum dot2i2iii (DQD) systems. From early theoretical work 
on the low-temperature properties of the two-impurity Kondo 
Hamiltoniani^ii it is known that localized moments form ei- 
ther each its own Kondo singlet with delocalized electrons in 
the leads - double Kondo phase (2K) - or they form local 
spin-singlet state (LSS) decoupled from delocalized electrons. 
The crossover between the two regimes is the consequence of 
competing energies of Kondo singlet with Kondo temperature 
Tjc and of LSS formation J. Similar results were obtained by 
the analysis of a two-impurity Anderson model by means of 
slave-boson formalismi^iA&Aii&IS and numerical renormaliza- 
tion group (NRG)22i2i. Resembling behavior was found also 
in particular regimes of a triple quantum dot system 2 ^. 

Here we focus on the role the inter-dot (capacitive) inter- 
action plays in the electron transport through serially coupled 
DQD. The influence of the capacitive coupling was already 
studied by means of the equation-of-motion method 23 ^ 2 ^ 6 ! 
which, however, fails to capture the Kondo correlations accu- 
rately. On the other hand, the Kondo correlations were consid- 
ered in the limit of vanishing tunnel coupling and strong inter- 
dot interaction where the ground state of the isolated DQD 
exhibits four-fold degeneracy: in addition to spin degener- 
acy also singly occupied states labeled with (1,0) and (0,1) 
are degenerate leading to orbital Kondo behavior 27 - 28,29,30 - 31 . 
Here (711,712) corresponds to occupancies n\ and 712 of the 
two dots. The simultaneous presence of spin and orbital de- 
generacy results in enhanced Kondo temperature Tk(N) ~ 
^/J'poN cxp(— 1/ J'pqN), where N is the degeneracy, /?o 
the noninteracting density of states, and J' the correspond- 
ing Schrieffer-Wolf prefactor 28 . The isospin Kondo resonance 
was observed indirectly in the measurements of enhanced 
conductance through DQDs2*i2iiiiA2i and carbon nanotubes^ 2 ^ 
and also directly, in the bulk, with techniques of the scan- 
ning tunneling microscopy 34 . Recently the enhanced Kondo 
temperature phase (EKT) was predicted for DQD with inter- 
dot interaction V = U in the absence of tunneling between 



the dots also at half-filling 



i.e., n 



ri2 = 2. The 



0.5 




3 

K 2 

1 





go 
in 



L -0.3 



-0.6 





/ \\ 7 \ 


„ ,-i 


I 
1 




\\ // 
\ \ ' / 




~ - N 7(=L0\ >! 


"(a)/ 


\A- - »/»/-H \-2t-J(v=o) \ \ 






is I >L 


■— v=u, t=r 








— V=U,t=TI5 








v=u, «=r/50 








— v=o, t=r 








-- v=o,t=r/5 










"(b) 


-i_ .,,,) J,. 








V 2K, EKT // 

s \ N. s* I - 












V 




\6ss/ 




e,U e,U 


"(c) 


\L£Sy 


'0 ' ' ' 'c 


-10 





10 


20 



(e+t//2+V)/T 



Figure 1: (Color online) Conductance (a), occupancy (b) and spin- 
spin correlation (c) of DQD for V = (thick), and V = U (thin). 
The vertical lines for V = U (dashed), and V = (dotted) corre- 
spond to transitions of isolated DQD between n = 0, n = 1 and 
n = l,n = 2 for t = T. Spin-singlet regimes (LSS, OSS) and 
Kondo regimes (2K, EKT, IK) are indicated with labels. Inset: DQD 
system with both: inter- and intra-dot Coulomb coupling, attached to 
noninteracting leads. 



ature Tjc) on top of an incoherent continuum in the density 
of states. However, the analysis of capacitively coupled DQD 
with finite inter-dot tunneling rates has been lacking to date 
and we address this issue in the present paper. 

We model DQD shown in inset of Fig. ^c) by the 
Anderson-type Hamiltonian H = Ha. + Hi, where Ha. cor- 
responds to the isolated dots 



Ha 



E 

»=1,2 



(etii + U UifTiii) + Vnin 2 - 1 ^(c| CT c 2(T + h.c. 



main characteristics of this phase at low temperatures is the 
enhanced width of the Kondo resonance (i.e., Kondo temper- 



with rii = riii + riii, rii a — c ia Cia-. The dots are coupled by a 
tunneling matrix element t and a capacitive V term. The on- 
site energies e and the Hubbard repulsion U are taken equal 
for both dots. H\ corresponds to the noninteracting left and 
right tight-binding lead and to the coupling of DQD to the 
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leads, 

Hi = -t 



i 5^0,1, 2 



cLci. 



where the sites are labeled as shown in the inset of Fig. Qc). 

In this paper we study the low temperature properties of 
DQD, derived from the ground state at T = 0, determined 
by the Gunnarsson and Schonhammer projection-operator 
method 36 ' 37 . The conductance is calculated using the sine 
formula 37 , G = G sm 2 [(E+ - E-)/4t L], where G = 
2e 2 /h and E± are the ground state energies of a large aux- 
iliary ring consisting of L non-interacting sites and an em- 
bedded DQD, with periodic and anti-periodic boundary con- 
ditions, respectively. The chemical potential is set in the mid- 
dle of the band, which corresponds to L electrons in the ring. 
The method proved to be very efficient for various systems 
with Kondo correlationsSiSi^iffi. Up to N ~ 5000 sites were 
used with variational ansatz based on ~ 100 projection oper- 
ators, which was sufficient to obtain converged results for the 
boundaries between various regimes. 

For large inter-dot tunneling t > U molecular bonding and 
anti-bonding orbitals with energies e b ( a ) = e^ft are formed. 
Whenever either of the two orbitals is singly occupied, the 
physics is analogous to the single impurity Anderson prob- 
lem with unitary transmission and Kondo correlations (IK 
regime). When the bonding orbital is doubly occupied the 
electrons form an orbital spin-singlet state (OSS) with dimin- 
ished conductance as the Kondo peak in the density of states 
is absent. The other limit, £ — ^ 0, is more interesting. When 
also inter-dot repulsion is moderate V < U the dots are ef- 
fectively decoupled and in the particle-hole symmetric point 
(e + U/2 + V = 0) the 2K state occurs. Regimes with dif- 
ferent occupancies are separated by the intermediate valence 
regime of width determined by the noninteracting hybridiza- 
tion r = t ,2 /t . 

Behavior of the system is less obvious in the intermediate 
region < t < U. In Fig. ^a-c) we plot the conductance, 
occupancy and spin-spin correlation for V = and V = U 
(T = to /25 and U/T = 15 is kept constant throughout the pa- 
per). The curves for t = V, V — U (dashed-dotted thin lines) 
are the most reminiscent of known results for single impurity 
case. The wide plateau in conductance curve corresponds to 
single occupancy regime. While the occupancy for the case 
with reduced t = T/5 (dashed thin line) does not change sig- 
nificantly, two maxima in conductance become discernible. 
One reflects the IK Kondo physics corresponding to the bond- 
ing orbital (n = 1 there). The other maximum (with unitary 
conductance) is known from the V — case 16 ' 21 and occurs 
generally at the transition between extended Kondo and de- 
coupled spin-singlet states at n ~ 2. For smaller t = T/50 
at n ~ 1 conductance is significantly reduced, but the ground 
state is still IK, which is signaled with the remaining plateau 
in occupancy and diminished fluctuations of the occupation of 
the relevant, (i.e., bonding or anti -bonding) orbital (not shown 
here). In this case the OSS state is not energetically favorable 
near n ~ 2, correspondingly the conductance peak with the 
unitary conductance is absent (likewise for V = 0, t = T/5, 
where the LSS is not energetically preferred). At n — 2 rather 
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Figure 2: (Color online) (a) Orbital Aiij = An„ and site Anf = 
Ari2 occupancy fluctuations at half-filling (n = 2) as a function 
of t/T for V = and V — U. Conductance (b) and spin-spin 
correlation (c), for V = (full) and V = U (dashed). 



the EKT phase is found, as discussed later. 

The V = case is seen to differ qualitatively from the V — 
U case. Crude insight into distinction between the two may 
be gained from treating just the filling properties of a detached 
system 40 . The first electron is added when e = t, and the 
secondwhene = -t+ J+ \{U+V)- \U-V\}/2, where J = 
[-\U-V\ + yf{U -V) 2 + 16t 2 } /2 is the difference between 
singlet and triplet energy of isolated DQD. When n = 2 the 
ground state is [a(\U) - JIT)) + /3(|20) - |02))]/V5, where 
a/ /3 = U/(V - U + y/ty - V) 2 + 16t 2 ). For V = and 
V = U case the values of e where the occupancy of a ground 
state changes from n = to n = 1 and from n = 1 to n = 2 
are indicated by vertical lines on Fig. Qa). The range of e 
where single occupation is favorable is diminished in the V = 
case. 

This observation remains valid also for DQD attached to 
the leads, Fig.^b). At first sight surprising feature is seen in 
Fig. [He), where the expectation value of a spin-spin correla- 
tor Si • S2 is seen to approach —3/8 near n — 2, but in fact 
it is just equal to the noninteracting result with two electrons 
occupying the bonding orbital. The fact that the spin-spin cor- 
relator approaches —3/4 near n = 2 in the V = 0, t = T case 
suggests that the local picture prevails. The electrons indeed 
form singlet in the local basis what is seen also in diminished 
value of the local occupancy fluctuations, Fig. |2ja). Addi- 
tional insight into precise role of inter-dot interaction can be 
obtained by re-writing the total Hamiltonian in the basis of 
orbital operators, 



Ha 



E 



u + v, 

£ a n a H — (n QT n Q x + n QT n a J 



+ 



+ V ^ n acrn blJ + ~ (Cflip - Sfljp) , 
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where notation a — b, b = a is used. The last term of 



and spin-flip 

(.sty 



consists of charge-flip Cff lp = T+T^ 
Sffip = S+S^ + h.c. operators, where 
are spin and Tr = c\-tC\\ — (T^)^ charge (isospin)il low- 
ering and raising operators for the orbitals A = b, a (or sites 
A = 1, 2). The full spin (isospin) algebra is closed with opera- 
tors S* = (n Q | —n a i)/2 andT^ = (n a — l)/2, respectively. 

When V = U, the spin- and isospin-fiip terms in are 
absent: the Hamiltonian is mapped exactly to the two-level 
Hamiltonian with intra- and inter-level interaction U with the 
bonding and anti-bonding levels coupled to even and odd 
transmission channels, respectively. When V ^ U this map- 
ping is no longer strictly valid. Taking V = case as an ex- 
treme example, we find two mechanisms responsible for this. 
Firstly, the electrons can avoid the inter-level repulsion by oc- 
cupying aligned spin-states in different orbitals, and secondly, 
charge flip-terms induce the fluctuations of charge between 
orbitals. Both mechanisms prohibit electrons to occupy the 
well-defined orbital states. 

In the rest of the paper we concentrate on the symmetric 
point of the model, n = 2. The manifestation of adequacy of 
orbital picture are the fluctuations of occupancy An^ = (n\ — 
(n\)) 2 , shown in Fig. |3a). For V = U orbital occupancy 
fluctuations An^(= An 2 for n = 2) are indeed smaller than 
local occupancy fluctuations An\ — An 2 , for all values of t. 
On contrary, for V = 0, An 2 < Aril onrv f° r * ^ 

To explain the conductance and spin-spin correlation shown 
in Fig.|2jb,c) we first note that the electrons bind into singlet 
(local or orbital, depending on the V/U), whenever J > J c = 
2.22V, which is known for V — results 13,43 , but we find 
that this result can be readily generalized for by taking the 
difference between singlet and triplet energies J also for V > 
0. This relation indicates that although the Kondo temperature 
is known to rise in the V = U case^ 5 . for n = 2, the OSS is 
energetically favorable for moderate t as J(V ~ U) oc t is 
larger when compared to J(V = 0) ~ At 2 /U . 

It should be noted that due to the variational nature of the 
method, the boundaries between various regimes are accu- 
rately reproduced. Sharp transitions of correlation functions 
in the vicinity of crossovers are less precisely determined. 
The transitions discussed here are actually smoothened into 
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Figure 3: (Color online) Phase diagram in (t/T,V/U) plane for n 
2. 




Figure 4: (Color online) (a) Fluctuations of charge in orbitals An„ = 
An b , (b) the spin-spin Si S2 and isospin-isospin correlation Ti -T2 
are plotted for t > t c (full), for t* < t < t c (dashed) and t < t* 
(dotted), (c) The conductance for the three cases. 



crossovers what happens generally when the parity break- 
ing terms (here the inter-dot tunneling) are present in the 
Hamiltonian 4 ™. 

For V — in absence of inter-dot tunneling t — > 0, each of 
the dots forms a Kondo singlet with its own lead (2K regime). 
When t is increased above some value t c , set by J > J c , 
the ground-state becomes a LSS, what is seen from increased 
value of —Si • S2 — > 3/4. The conductance is small in both 
cases due to the effectively decoupled leads in the 2K case 
or deficiency of states near the chemical potential in the LSS 
case. The conductance G — > Go at some point what is most 
vividly understood in terms of the phase shift between odd and 
even transmission channels, changing from to ir during the 
transition 16 . For V < U the behavior does not change qualita- 
tively, as the transition remains close to 2Tk ~ J- In the case 
V ~ U numerical results show deviations from this predic- 
tion: the transition occurs at smaller values of t than set with 
relation J > J c . We explain this by increased Kondo tem- 
perature in the EKT phase as suggested from NRG results 35 . 
Here the transition between 'SU(4) Kondo' 44 and OSS occurs 



St/(4) 
K 



J(V = U) = 2t* c 



atT, 

Further increasing V above U we find, that in contrast to 
the strict t = case where the ground state for V > U is a 
non-Fermi liquid charge-degenerate phase 35 , a finite t breaks 
the degeneracy resulting in a (Fermi liquid) phase with in- 
creased magnitude of isospin-isospin correlation Ti • T2: a 
local charge singlet (LCS). LCS corresponds to the ground 
state of a detached system for (V — U)/t 3> 1. 

In Fig. |3 the phase diagram of DQD for n = 2 is shown. 
The border line between Kondo and local-singlet phases (full 
line with diamonds) is characterized by unitary conductance 
and the abrupt increase of —Si • S2 for the transition between 
2K and LSS and of — Ti • T2 for transition between EKT and 
LCS [see also Fig. |2](a,b) and Fig. Etb,c)]. In the t, t' 
limit the states (2, 0), (0, 2) and four states (1, 1) are degener- 
ate, and (— Ti ■ T2) = 1/6, if the expectation value is taken 



4 



with respect to the ground state where each of these states 
occurs with equal probability. We quantitatively characterize 
the transition to the EKT phase (full line with circles) when 
— Ti • T2 > 1/8. The border lines which we use to distin- 
guish between local and orbital regimes (dotted) are given by 
An 2 b = Anl for V < U and by -Si • S 2 = 3/16 (half of 
noninteracting value) for V > U. 

In Fig.|4]we sample the phase diagram at different t by in- 
creasing V. We plot three different cases. In the first case 
(full curve) DQD is in LSS already at V = 0. By increasing 
V we see a crossover to LCS through the orbital regime. In 
the other two cases we start in 2K. Now if t < t* (dotted) we 
stay in the Kondo regime until the transition to LCS happens 
with the unitary conductance at the transition point. However 
if t > t c (dashed) there is an area in between where LSS is 
energetically favorable. From the point where the system is in 
LSS, it follows the same scenario as in the first case. 

Finally, let us briefly comment on V > U part of the phase- 
diagram. This regime could arise when the dots are capac- 
itively coupled while the on-site repulsion is reduced due to 



the coupling to local Einstein vibrational modes**ai. Here the 
competition between LCS and EKT is relevant for the ground- 
state. Our numerical data confirms NRG predictions that the 
Kondo regime is relevant only for small increase of V over U. 
For parameters used here, we find critical V c ~ U + [7/20. 

In conclusion, we have presented the phase diagram of a 
pair of tunneling-coupled quantum dots with Coulomb inter- 
and intra-dot interaction. In general, the behavior of V ~ U 
case was found to differ from the V ~ case due to differ- 
ent arrangement of orbital levels of a decoupled system. In 
addition, for V ~ U the correct exchange coupling is not 
J ~ At 2 /U but rather J ~ 2t, hence the spin-singlet bind- 
ing energy is increased. The increase in J is seen to sur- 
mount the enhancement of Kondo temperature which is due 
to the additional degeneracy in the charge sector. Therefore 
this enhanced Kondo effect for n = 2 should be discernible 
in measurements only in a tiny window of tunneling rates and 
inter-dot vs. intra-dot interaction. 

We thank R. Zitko for valuable discussions and acknowl- 
edge the support of MSZS under grant Pl-0044. 
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